The redshift distribution of gamma-ray bursts collected in 
the BATSE Catalog is compared with the star formation rate. 
We aim to clarify the accordance between them. We also study 
the case of comoving number density of bursts monotonously 
increasing up to redshift ^ (6 - 20). A method independent of 
the models of the gamma-ray bursts is used. The short and the 
long subgroups are studied separately. The redshift distribution 
of the long bursts may be proportional to the star formation 
rate. For the short bursts this can also happen, but the propor- 
tionality is less evident. For the long bursts the monotonously 
increasing scenario is also less probable but still can occur. For 
the short bursts this alternative seems to be excluded. 
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Abstract. 



1. Introduction 

The Burst and Transient Source Experiment (BATSE) 
instrument on the Compton Gamma Ray Observatory 
( |Meegan et al. 200' ll detected around 2700 gamma-ray bursts 
(GRBs). From this data set it follows that there are two different 
subgroups of GRBs, "short" and "long" ones (see, e.g., Balazs 
et al. (2003) and the references therein); the existence of fur- 
ther subgroups is not excluded (Horvath 1998 Horvath 2002 
IHakkila et al. 2003l|Eorgonovo 20 04). 

The redshift (z) distributions of the GRB subgroups are not 
satisfactorily known. For the long subgroup the redshifts are 
known for ~ 40 long GRBs from the afteglow measurements 
provided after 1997. On the other hand, for the short GRBs only 
one direct redshift is known (Greiner 2005 1. In the case of the 
further subgroups even the physical reality of these subgroups 
is in doubt ( Hakkila et al. 2003d . (For more details about these 
questions see, e.g., the surveys of Meszaros (2001) and Piran 
(2004).) 

Here we concentrate on the question: Can the redshifts of 
GRBs be distributed in accordance with the distribution of the 
other objects arising in the star formation regions? 

This question is discussed, e.g., by Porciani & Madau 
(2001), and the question is answered in essence positively. 
Nevertheless, great care is still needed for three reasons. First, 
Porciani & Madau (2001) note that - assuming different lu- 
minosity functions for GRBs - one can always be consistent 
with the observations. Thus, a study having minimal theo- 
retical assumptions about the luminosity functions would be 
useful. Second, once one is in accordance with the star for- 
mation rate (hereafter SFR), an alternative distribution can 
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also fit the data. For example, several papers (Meszaros & 
Meszaros 1995, Meszaros & Meszaros 1996, Horvath et al. 
1996, Reichart & Meszaros 1997, Hogg & Fruchter 1999, 
Norris et al. 2000, Schaefer et al. 2001, Lloyd-Ronning et al. 
2002, Non-is 2002, Bago ly et al. 2003|lLin et al. 2 00T) suggest 
that, alternatively, a permanent monotonous increase of the co- 
moving number density of GRBs up to z - (10 - 20) is also 
possible. For example, Hogg & Fruchter (1999) consider this 
scenario as an acceptable alternative. Hence, this alternative 
should also be studied. (Note that Hogg & Fruchter (1999) 
consider also a third possibility: The density of GRBs is pro- 
portional to the total integrated stellar density. This scenario is 
found to be less probable, and therefore it will not be discussed 
here.) Third, the evidence has increased for the intrinsic diffe- 
rence between short and long GRBs (see Balazs et al. (2003) 
and references therein). But, it is not sure that both subgroups 
are in accordance with the SFR. Hence, a new comparison of 
GRB redshift distributions and star formation scenarios sepa- 
rately for the short and long GRBs is needed. In this paper a 
procedure will be provided for the two subgroups separately 
with minimal theoretical assumptions about both the SFR sce- 
nario and its alternative. 

The paper is organized as follows. In Sect. 2 we discuss the 
method. Sect. 3 contains the theoretical calculations, compared 
with the observations from the BATSE Catalog. Sect. 4 veri- 
fies these theoretical calculations. Sect. 5 discusses the obser- 
vational biases. Sect. 6 gives the comparison of theory with ob- 
servations. In Sect. 7 the results are summarized. 

2. The method 

Balazs et al. (2003) studied the fluences and the durations of 
GRBs, collected in the BATSE Catalog ( |Meegan et al. 2001| l. 
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For the fluence F it is fulfilled 



F = 



= c{z)Eis, 



(1) 



where Eiso is the total emitted energy of the object assuming 
isotropic emission, and where diiz) is the luminosity distance 
of the object at redshift z- Balazs et al. (2003) show that log F 
is distributed normally separately for the two subgroups, and if 
c{z) and E, ™ are independent variables, then one has 



(2) 



where cr is the dispersion for the given quantity denoted by the 
index. The independence seems to occur both for the short and 
long subgroups, respectively. The statistical fitting of log F al- 
lows - separately for the two subgroups - a normal distribution 
with a dispersion o-\ogF^ and if o"iogf(j) and criog£„,,, are assumed 
to be comparable, then the Cramer theorem says that both vari- 
ables log Eiso and log c(z) should also be distributed normally - 
again separately for the two subgroups. It is also possible that 
either crxogdz) » o'iog£,„ or criogcfe) <k o"iog£,„, occurs. If this 
is the case, then nothing can be said about the distribution of 
the variable having a much smaller dispersions - it may be nor- 
mal; however, the variable with dominating dispersion must be 
normal. In Balazs et al. (2003) it is argued that the condition 
Ciogcfe) » o"iog£;„ is unlikely. Hence, only the normal distrubu- 
tion for log Eiso is a reasonable conclusion. No detailed discus- 
sion about the distribution of c(z), i.e. about the redshift distri- 
bution of GRBs, is provided by Balazs et al. (2003). 

In this article we will concentrate on the redshift distribu- 
tions of the two subgroups. 

The key idea of this article is the following. 

We assume that a given subclass of GRBs in the BATSE 
Catalog is distributed in accordance with the redshift distribu- 
tion of the objects in the star formation regions. We then calcu- 
late the theoretically expected distribution of log c(z), and the 
theoretically expected dispersion o-\ogc(z),theor- These theoretical 
calculations are general cosmological ones, and hold for any 
subsample. The theoretical dispersion can be compared with 
the observed crxogp- Because we have two different observa- 
tional values of cr\o^f for the two subgroups, we can well make 
this compare the two subgroups separately. We will use one 
o"iog cfe),/;icor value twice - for the two subsamples of GRBs. 

It is a necessary condition that o-\ogc(zuheor < o"iogF be ful- 
filled. If this condition is not fulfilled, then the primary assump- 
tion must be rejected. For a given subgroup the condition may 
be fulfilled, but not fulfilled for the second one, because for the 
second subgroup the observed crjog p values are smaller. 

If o"iog c(z),//if or ^ cTiogf holds, then the primary assumption 
is acceptable, because in this case log c(z) can have an arbitrary 
distribution. If <J\ogc{z) and cr\ogE,„, are comparable then both 
variables should be distributed normally. For log Ejso this is true 
( IBalazs et al. 200 3 1. But for logc(z) the requirement for nor- 
mal distribution is artificial. Nevertheless, it can happen that the 
distribution "mimics" a normal distribution. This means that - 
mathematically - the functional behavior of the distribution of 
log c{z) is similar to a Gaussian curve. If this nearly Gaussian 
distribution of log ciz) occurs, then this function together with 
the normal distribution of log Eiso allows the acceptance of the 



primary assumption. If criog^,^ <k cr\ogc(z),tiwor < o"iogF holds 
then the situation is similar; the primary assumption is accept- 
able, when the distribution of log ciz) mimics a Gaussian curve. 
This procedure does not need any assumption about the model 
of the GRB. 

If one assumes an accordance of SFR and the occurence of 
GRBs, one has to specify a SFR for the calculation of the dis- 
tribution of c(z). The SFR is taken from the literature. Several 
papers (Madau 1995 Mada uetal. 19961 ISteidel et al. 19991 
IPorciani & Madau 2001 Wilson et al. (2002), Tonry et al. 
(2003), Giavalisco et al. 2004a, Giavahsco et al. 2004b, Dahlen 
et al. 2004, Strolger et al. 2004) suggest the following behavior 
of the comoving number density for the star forming regions: 
1 . At z ^ 1 the density should be greater than at z ^ 0. 1 roughly 
by ^ 10 times; 2. At higher redshifts, up to z ^ 6, the comoving 
number density should remain nearly constant, or there should 
be a weak decrease. Therefore, there should be four indepen- 
dent parameters in the description of the SFR: the two typical 
redshifts (zbreak, where the density is peaked; Zma.t, up to which 
the rate is defined), and the two characteristic density ratios 
([(density at the peak redshift)/(density atz = 0)] and [(density 
at ZmoA )/(density at z = 0)]). There is also a fifth parameter, the 
density at z = 0, but this value will always disappear from our 
calculations. The exact functional form of the SFR may be ana- 
lytically expressed by different empirical functions. For exam- 
ple, in Strolger et al. (2004) a smooth four-parametric function 
is used. Nevertheless, other functional dependences may also 
be chosen that contain the four free parameters. We will take - 
in accordance with Wilson et al. (2002) - a power-law depen- 
dence of the form oc (1 h- z)°' between z = and Zbreak, and 
then oc (1 -H z)^- between Zbreak and Zma.r- Reasonable values of 
Di are roughly between 2 and 4, and for D2 between - 1 and 0. 
The reasonable values of parameters Zbreak are around 1.0-1.5 
and Zmax around 4-6. These values define the range of param- 
eters in the theoretical calculations for o-\ogc(z).iiieoi-- Of course, 
all reasonable cosmological models should also be considered. 

No "identical" number density is needed because the nor- 
malization constant - the fifth parameter - vanishes. Thus, under 
the "accordance", we simply mean a "proportionality". 

We will also discuss the alternative scenario. For this al- 
ternative we take in the theoretical calculations a monotonous 
growth of the comoving number density of GRBs up to very 
high redshifts. Mathematically, this is obtained from the SFR 
density, if in it Zmax = Zhreak is taken. Then we have a oc (1 -nz)^' 
behavior between z = and z = Zmax, where Zmax - (5 - 20). 
The parameter D2 does not exist in this case. Thus, we have 
only two free parameters: Zmax - Zbreak and Di. Having one 
single theoretical value, this value should be compared twice 
with the observed criogF from the Balazs et al. (2003) values 
separately for the two suub groups. 

Let the real physical density of GRBs at the redshift z be 
given by n(z). Its unit is Mpc 'yr ' . Then the number of GRBs, 
being at the infinitesimal redshift interval [z, (z + dz)] and ob- 
served by an observer at z = 0, is given by 



N(z) dz 
= n{z) 



Andljz,Ho,nM,n_A)_ ^ -cdt{z,Ho,nM,nA) )j^^ 



(l+zV 



(3) 
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where c is the velocity of light, H„ is the Hubble con- 
stant, Qm is the ratio of the density of matter to the criti- 
cal density, Qa is the dimensionless cosmological constant. 
The function dt{z,H„,D.M,0.i^ldz is dependent on z and on 
the three cosmological parameters (//„, Q.m, ^\)- The distance 

dpM{z,Ho,aM,^K) = (c///„)e(z,QM,nA) = ddz)l{\+z) is 

the proper-motion distance. This formula is a standard cosmo- 
logical one (Carroll et al. 1992). 

Mathematically, N{z = 0) = 0, because dpMiz = 0) = 0. In 
addition, dpM remains finite, if z ^ oo. Hence, to have N{z) - 
at the limit z ^ °° one needs a slower than n{z) oc { \ + z)^-^ 
increase for z ^ oo. In other words, if n(z) remains finite, or 
runs to oo for z — > oo but slower than oc (1 -t- z)^'^, then the zero 
limit at infinite redshift is ensured for N{z)- This is always the 
case here. This global behavior of N{z) is independent of the 
cosmological parameters, and hence may "mimic" a Gaussian 
distribution both for the SFR and for the alternative scenario. 

n{z) is the real physical density of GRBs under the con- 
dition of isotropic emission. Nevertheless, this assumption is 
not a loss of generality. Assume that GRBs emit in the solid 
angle cj (in steradian), where o) < An. Then the real phys- 
ical density of GRBs is given by n{z){A-nl<jS), but the ob- 
served number of GRBs is lowered by the factor w/(47r), and 
hence all formulae about the observed quantities hold. Thus, 
our theoretical calculations are independent of the size of 
the beaming JFrail et al. 20 01 Amat i et al. 2"003l IPiran 20041 
IGhirlandaetal. 2004t . 

Having N{z), it is straightforward for any function f{z) 
to define its mean f{z\,Z2) and dispersion cr^ in the interval 

[Z\,Z2\ (Zl < Zl). 

A widely accepted cosmological model is given hy Q.m - 
0.3, Qa = 0.7, where Q.m + Ha - 1 holds exactly 
( |Tonry et al. 2 003 1. Nevertheless, from the observational point 
of view, the case Qa - with Qm - (0.2 - 0.4) cannot be ex- 
cluded (see Meszaros (2002), Tonry et al. (2003) and references 
therein). Therefore, for maximal generality, we should discuss 
both cases. We will also discuss the case Qm = 1 with Qa = 0. 
This eventuahty is rejected by observations, but it can serve as 
a verificator. In the case with Qm + Qa = 1, if Qa is decreased 
from 0.7 to zero, one should obtain this limit. In the case with 
Qa = the value Qm can be increased from 0.2 to 1, and one 
should again obtain this limit. Hence this case can be taken as 
the limit for verificator. We proceed similarly to Meszaros & 
Meszaros (1995), Meszaros & Meszaros (1996) and Horvath 
et al. (1996), where this simplest model was also discussed. 
The cosmological parameter Hg need not be specified for our 
calculations. 

3. Expected values for (T\ogc{z),theor 

Now assume that n(z) is in accordance with SFR. Then n(z) 
may be described by 

n{z) = no(l -H Z)^(l + Zf' , < Z < Ztreak , 

n(z) = nod + Z)^^^-(1 + Zbreakf''"-, Zbreak <Z< Zmax , (4) 

where Zbreak ^ (1.0 - 1.5), z,„„ ^ (4.0 - 6.0), ^ (2 - 4) and 
D2 - (-1 - 0), and where «(z) = for z > Zmax- 



The extra (1 + z)^ factor in the definition of n(z) is needed 
here because in our formulae «(z) means the real physical 
proper density, not the comoving one. In the general case, one 
must calculate the integrals numerically. The final cr\ogc(:),theor 
will depend on the four parameters: Zbreak, Zmax, Di and D2 . The 
value Ho vanishes from the formulae, and need not be specified. 

This behavior may easily be changed into the alternative 
scenario with a monotonous increase of the comoving den- 
sity with Zhreak - Zmax (i-B. no sccond part with the exponent 
D2), and for Zbreak - Zmax One can take any value between 
z = 5 and z — 20. Then the best value for Di seems to be 
0.5 ( Reichart & Meszaros 19971 1 with a small scatter between 
Di = OandDi = 1. 

Using the equations from the previous Section, one has for 
Qm = 1 and Qa - 

c(z) = [An {clHof (1 + z) Qizfr' = [d + z) Andl^V , (5) 

where Q(z) - 2(1 - (1 -H z)"'^^). Hence, calculating the dis- 
persion of log c(z), it is enough to calculate the dispersion of 
log[(l + z)Q{zflA] = log( VTTI - 1)2 = 2 log( Vm - 1). 

The results for o-\ogc(z),theor. with several characteristic Zbreak, 
Zmax, Di and D2 values are shown in Table 1. 

o"iog c(z),//icor. does not depend strongly on the parameters 
Zbreak, Zmax, D\ and D2, becausc all values lead to the values 
(0.41 - 0.58). 

The change given by the different Zmax is negligible. In 
addition, if we change Zmax to a larger interval (say, between 
2 < Zmax < 10; not written in the Table 1), the changes remain 
small. This is because N{z) is decreasing at z > 2. 

Similarly, the changes given by the different values of Zbreak 
are also small. If we change its value at a larger interval (0.75 < 
Zbreak < 2.0; not Written in Table 1), the changes remain small. 
This is also expected from the behavior of N{z), because it in- 
creases up to z - 1, and then it decreases. The change of Zbreak 
influences this behavior only weakly around the maximum of 
N{z). 

There is a clear tendency for D\ and D2: increasing D\ de- 
creases the dispersion, and, similarly, decreasing D2 also de- 
creases a-\ogc(z).iheor.- Greater D\ and/or smaller D2 strengthens 
the "peakness" of N{z), and thus decreases the dispersion be- 
cause more objects are concentrated around ztreak, where the 
mean should also be. Thus, one expects smaller dispersion. 

We also calculated the values for the alternative scenario, 
when there is a monotonous n(z) oc (1 -1- z)^' increasing up to 
very high redshift Zmax- Similar dispersions were obtained with 
a very weak dependence on Zmax and D\ . 

For the case < Qm < 1, Qa = one may use 
Q{z, Qm) = 2 [2 + Qm(z - 1 ) - (2 - Qm) VI+^mz)] / [^^i l+z)]. 
The key difference, compared to Qm - 1, follows from the 
smaller value involved here. We also probed several different 
< Qm < 1 (not given in Table 1), and obtained a clear in- 
crease of cr\ogc(z),iheor with Qm. Otherwise the trends are iden- 
tical to that of Qm = 1 . The alternative scenario shows similar 
trends with larger values. 

For Qm = 0.3, Qa = 0.7 we need to determine the values 
of dpM numerically. The results are similar to that of Qm - 
0.3, Qa - 0. Thus, the values for Qm = 0.3 are similar in 
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Table 1. The expected dispersion criogc(z),th£or for the different 
parameters in n(z). 



Dl D2 Zbreak Z„, 



— 1 
Qa = 



Q.M - 



0.3 

= 



Qm = 0.3 
JIa = 0.7 













'-^log c{z)-lheoi'. 




2.5 


0.0 


1.0 


4.0 


0.51 


0.35 


0.39 


2.5 


0.0 


1.0 


6.0 


0.58 


0.40 


0.40 


2.5 


0.0 


1.5 


4.0 


0.46 


0.31 


0.36 


2.5 


0.0 


1.5 


6.0 


0.53 


0.35 


0.38 


2.5 


-1.0 


1.0 


4.0 


0.51 


0.35 


0.42 


2.5 


-1.0 


1.0 


6.0 


0.58 


0.42 


0.42 


2.5 


-1.0 


1.5 


4.0 


0.46 


0.32 


0.38 


2.5 


-1.0 


1.5 


6.0 


0.52 


0.37 


0.39 


3.0 


0.0 


1.0 


4.0 


0.49 


0.34 


0.39 


3.0 


0.0 


1.0 


6.0 


0.56 


0.40 


0.41 


3.0 


0.0 


1.5 


4.0 


0.44 


0.30 


0.35 


3.0 


0.0 


1.5 


6.0 


0.51 


0.35 


0.38 


3.0 


-1.0 


1.0 


4.0 


0.51 


0.36 


0.42 


3.0 


-1.0 


1.0 


6.0 


0.54 


0.41 


0.43 


3.0 


-1.0 


1.5 


4.0 


0.44 


0.31 


0.37 


3.0 


-1.0 


1.5 


6.0 


0.50 


0.34 


0.39 


3.5 


0.0 


1.0 


4.0 


0.48 


0.33 


0.39 


3.5 


0.0 


1.0 


6.0 


0.55 


0.39 


0.42 


3.5 


0.0 


1.5 


4.0 


0.42 


0.28 


0.34 


3.5 


0.0 


1.5 


6.0 


0.49 


0.34 


0.38 


3.5 


-1.0 


1.0 


4.0 


0.48 


0.35 


0.42 


3.5 


-1.0 


1.0 


6.0 


0.55 


0.40 


0.44 


3.5 


-1.0 


1.5 


4.0 


0.41 


0.29 


0.35 


3.5 


-1.0 


1.5 


6.0 


0.48 


0.35 


0.39 


0.0 




6.0 


6.0 


0.48 


0.53 


0.31 


0.0 




10.0 


10.0 


0.54 


0.58 


0.33 


0.0 




20.0 


20.0 


0.59 


0.63 


0.34 


0.5 




6.0 


6.0 


0.47 


0.49 


0.33 


0.5 




10.0 


10.0 


0.49 


0.54 


0.36 


0.5 




20.0 


20.0 


0.51 


0.59 


0.43 


1.0 




6.0 


6.0 


0.47 


0.45 


0.34 


1.0 




10.0 


10.0 


0.52 


0.48 


0.40 


1.0 




20.0 


20.0 


0.56 


0.52 


0.50 



both cases, and depend weakly on the value of the cosmological 
constant. For the alternative scenario we obtain smaller values. 



4. K-correction at gamma range 

Before comparing the theoretical values of Table 1 with the 
fit of Balazs et al. (2003) we still discuss the values of Table 
1. This is needed because in the previous considerations the 
assumption of the independence of c(z) and Ejso was crucial. 
Because c(z) is a function of z, this assumption is equivalent to 
the assumption that Ejs„ does not depend on z - a non-trivial 
assumption. 

In our case Eiso means the total emitted energy of the pho- 
tons having energies Ei{l + z) < E < £2(1 + z) under the 
assumption of isotropic emission. But this Eis„ can be different 
to that of the remitted total energy of the GRB in the interval 
Ei(l + z) < E < £2(1 + z) because we measured the fluence 
and the spectrum at the interval E\ < E < E2, and from this 



measurement we made extrapolations about the emitted energy 
at the interval shifted by a factor (1 + z). Hence, Ei^o used in 
Eq.l may depend on z even in the case when the real bolomet- 
ric isotropic energies - emitted by GRBs - do not depend on z. 
In essence, we have the analogy of K-correction in the gamma 
range. 

This question is discussed by Bloom (2003) for GRBs, and 
here we will proceed similarly. 

The case Eiso = Eisoiz), where Eiso is the quantity appearing 
in Eq.l, may be reformulated as 

Eiso(z) = Eiso(0)qiz) , (6) 

where E/soiO) is the emitted energy at z = between the pho- 
ton energies E[ and E2, and q(z) is some function of z to be 
determined from the time integrated spectra of GRBs. Thus, 

(l+z)Eiso(z) 



F = 



: c(z)qiz)EisoiO) . 



(7) 



47Tddz)^ 

This means that one should simply consider q(z)c(z) instead of 
c(z) itself, because then q{z)c(z) and EisoiO) are independent 
variables. The effect of K-correction is included into the func- 
tion c(z), and - to be the most accurate - one has to estimate the 
theoretically expected values of c^^ag^g^z)o(z)), instead of the val- 
ues CTi^gciz)- We estimate from the spectra of GRBs how much 
the value of q(z) can be changed between z = and, say, z = 20. 
From its definition it follows that, mathematically, q(z) is a di- 
mensionless non-negative function of z, and q(0) = 1 . 

Because we consider the whole emitted energy during the 
existence of a GRB, we should consider the time averaged 
spectra of GRBs. This can be represented by the semi-empirical 
Band spectrum JBand et al. 1993t taking the form 

S(E)^Ai(E/lOOTexp(-E/Eo), E<(a-/3)Eo, (8) 

S(E)^A2(E/100f, E>{a-/3)Eo, (9) 

where A2 takes the values ensuring the continuity at 
E = (a - P)Eo, and where E and are in keV. Then 
S (E)dE defines the number of received photons accross crn^ 
and having energies in the interval [E,(E + dE)]. Several 
papers (|Ryde & Svensson 2000| |Ryde & Svensson 2002| 
P^yde & Petrosian 2002 1 suggest that this form decribes the 
known GRB spectra; only the three parameters a,j3, Eo should 
be varied. For long bursts the typical values are a = -1, 
yS = -2 and E^ = 150 feV (Band et al. 1993 1. It is also known 
that this behavior can be extrapolated up to photon energies 
of tens of MeV dTavani 1996 1. Remarkably, this behavior of 
spectra holds also for short bursts (Ghirlanda et al. 2004 1; the 
only difference is given by the fact that the oc E" exp{-E/Eo) 
behavior holds in the discussed range leading to a much harder 
spectrum; hence Eo should be much higher than hundreds of 
keV, and the steeper slope with oc E^^ should not occur. The 
essentially different behavior of the short and long GRBs at the 
range > 300 keV was recognized also by Balazs et al. (2004). 

Taking o- = -1 and /3 - -2, and taking into account that 
for energies one should take ES (E), one obtains an analytical 
formula for q(z) in the form 



q(z)^ 



■(l+z)& 



S(E)EdE 



VS{E)EdE 
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_ exp(-xi(l +z))-e ^ + e ' In X2 + e ' ln(l + z) ^^^^ 
exp(-xi) - + e~^lnx2 

where x\ = Ei/Eg, and X2 = E2/E0. This formula holds for 
£■1(1+2) < Eo, if Eo < E2- For higher z > Zo one obtains 
q{z) - const - q(Zo), where I + Zo - Eg/Ei. This formula 
can be used to estimate the behavior of qiz) for the long bursts, 
because Eg < E2 is fulfilled for long GRBs. On the other hand, 
for short burst one has E,, > E2, and hence one may write 

c:s:s(E)EdE 

J,;S(E)EdE 

^ exp(-xi(l+z))-exp(-x.(l+z)) ^ ^ ^^^^ 

exp(-xi) - exp(-X2) 

because exp(-x2) <^ exp(-xi), and exp(-x2(l + z)) 
exp(-xi(l + z)) hold. 

In the case of short GRBs one has xi < 25/1000 = 0.025, 
and hence up to z = 4, q(z) is between 1.0 and 0.9; i.e. the 
change is smaller than 10%. For long bursts, substituting the 
different values of z with x, = 1/6 = 0.167, X2 = 1000/150 = 
6.67 one may verify that up to z - 5 the change of qiz) is 
smaller than 16%. Hence, assuming the accordance of Niz) 
with SFR, the predicted dispersions in the previous Section 
are not influenced more than (10-16)%. The effect should be 
smaller for short bursts. Because the parameters a,^, Eg may 
vary, more uncertainty should also be allowed, but even this 
variation in the parameters does not change the situation essen- 
tially, because the variation of parameters - on average - should 
cancel. Therefore, the estimation obtained for the typical val- 
ues of a,l5,Eo should hold. Hence a maximum ^ (10 - 16)% 
uncertainty can be expected from this effect; it wiU be smaller 
for short GRBs. Bloom (2003) in an other discussion of this 
effect also allows a ^ 20% uncertainty. All this means that the 
values of Table 1 with Zbreak - 1 can be used, but one should 
keep in mind that a 20% uncertainty may occur 

If one assumes «(z) oc (1 -h z)^' up to the very high red- 
shifts, then for the long GRBs the situation will be even better, 
because for the very large redshifts (up to z = 20 or so) the 
effects of K-correction should be even more negligible due to 
the behavior S (E) oc E~^ at the range of a few MeV keeping 
q{z) constant. Contrary to this, for short GRBs this behavior 
need not be fulfilled, because for them there is no indication 
for 5(£) oc £-2 around IMeV from Ghirlanda et al. (2004). In 
other words, if the short bursts are at z - (5 - 15) in accordance 
with n(z) oc (1 -H z)^', then for them the K-correction will not 
be negUgible, because q{z) will differ from one (it should go to 
zero as z increases). 

5. Biases 

In Sections 2-3 we have calculated the theoretically expected 
dispersions cr\o^c{z-),theon and in the previous Section we argued 
that the values of Table 1 - with some care in the case of the 
alternative monotonous growth for short GRBs - can be used as 
the expected theoretical dispersions in the comparison with the 
observed cxiogir. But it is well-known that the fluence itself is a 
biased quantity in the BATSE Catalog (see Balazs et al. (2003) 



for more details). Hence, the values of Ei^o can be biased. The 
theoretical values of criog c(z),theor from the previous Sections can 
be used, but the observed values of criogF must be taken with 
care. 

There are two different types of biases. The first type is 
given by the BATSE threshold leading to the fact that some 
faint GRBs are simply not detected. Then it is - in principle - 
possible that if they were also detected, criogF would be differ- 
ent. The second type of biasing follows from the fact that the 
fiuences themselves are affected by systematic uncertainities, 
and if they were free of them, again it would be possible that 
ciogF would be different. 

These biases of the fluence are discussed in detail by Balazs 
et al. (2003). It is shown there that the first type may be over- 
came quite easily: it is enough to study only the bright GRBs. 
For example, if only the bursts with > 1 photon/(cm2s) are 
taken {P^ is the peak-flux on the 64 ms trigger), then this trun- 
cated sample of bright GRBs should not be influenced by the 
threshold effects, and hence also by the first type of bias. 

A problem can arise from the fact that this truncated bright 
sample might not represent the whole BATSE sample. If the 
value of cTiogf is used from the whole sample, then this value 
may not be precise due to the bias; if the value of criogf is used 
from the truncated bright sample, then this value is bias-free, 
but it may not represent the value of the whole sample. 

This controversy in the general case would need a detailed 
study. Furtonately, in some special cases it is simply solvable. 
For example, this does not exist in the special case when the 
obtained values of o"iogF from both samples are identical. 

For the whole sample of long GRBs, as it is given in Balazs 
et al. (2003), one has (XiogF - 0.66. For the truncated bright 
subsample with £"64 > 1 photon/(cm^s) we obtained from a new 
calculation criog^ - 0.66. Also for the short GRBs we obtained 
from the new calculations the same crxo^p = 0.58 value both 
for the full sample and for the truncated sample with Pm > 
1 photon/(cm2s). Hence, the special case of indentical values 
is usable here. To check these identical values we calculated 
o"iogF for several limiting £"64- We obtained the behavior that 
criog/7 does not change for the truncated samples with P^ > 
(0.2-2.5) photon/(cm2s); for even bigger P^ there is a moderate 
decrease. This holds for both subgroups. Thus, the bias of the 
first type is unimportant in this article. 

The second type of bias was shown to be unimportant by 
Balazs et al. (2003). 

6. Comparison of theoretical predictions with 
BATSE fluence dispersions 

For the long GRBs, if one assumes that they are distributed 
in accordance with SFR, it follows that cr\o^c(z),theor = (0.28 - 
0.42) < 0.66 for = 0.3 independently of the value of the 
cosmological constant, and practically independent of the four 
parameters of SFR. Taking into account the uncertainty from 
the K-correction (^ (16 - 20)%) the condition a\o^c(z),theor < 
(TiogF holds. Because, in addition, the distribution of log c(z) 
can mimic a normal distribution, the redshift distribution of 
long GRBs can be in accordance with the SFR. 
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For the short GRBs, the requirement of o-iagc(,xtheor — 
(0.28 - 0.42) < 0.58 holds, but it is less obvious, because of 
the smaller values of the measured o-j^oF- The occurence of 
o'iog c(z),iheor - CTiogF cannot be excluded for Qm - 0.3 (in- 
dependently of the value of the cosmological constant). The 
relatively smaller K-correction (^ 10%) does not change the 
situation. All this means that, contrary to the long GRBs, here 
the accordance with SFR is less certain. 

If one considers the alternative scenarios, then 
Clog c(z),theor - CTjog /? cannot be excluded for = 0.3 
for the long GRBs with zero cosmological constant (with 
non-zero cosmological constant the situation is better). On 
the other hand, if one considers Di = 1 and takes the whole 
BATSE sample with (TiogF - 0.66, then one can still arrive 
in an accordance. The situation was less satisfactory than for 
long ones assuming accordance with SFR. 

For the short GRBs the validity of the alternative sce- 
nario is practically excluded because for any reasonable case 
o'\ogc(z),iheor - ClogF - 0.58 Can occur In addition, here also 
the K-correction may not be negligible. 

7. Conclusions 

Comparison of o-\o„c(z),theor with the observational criogf gives 
the following results: 

1 . The comoving number density of the long GRBs in the 
BATSE Catalog may be proportional to SFR. 

2. The comoving number density of the short GRBs in the 
BATSE Catalog can still be proportional to SFR, but the situa- 
tion is less obvious than in the case of long GRBs. 

3. For long GRBs in the BATSE Catalog the monotonous 
increase of the comoving number density in the form oc ( 1 +z)^' 
up to z ^ (6 - 20) can still occur, but is less clear than the case 
when this density is proportional to SFR. 

4. For the short GRBs in the BATSE Catalog, similar 
monotonous increase is excluded. 
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